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Abstract 



A model describing the evolution of a liquid crystal substance in the nematic 
phase is investigated in terms of two basic state variables: the velocity field u 
and the director field d, representing the preferred orientation of molecules in 
a neighborhood of any point in a reference domain. After recalling a known 
existence result, we investigate the long-time behavior of weak solutions. In 
particular, we show that any solution trajectory admits a non-empty w-limit set 
containing only stationary solutions. Moreover, we give a number of sufficient 
conditions in order that the w-limit set contains a single point. Our approach 
improves and generalizes existing results on the same problem. 
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1 Introduction 

In this paper we analyze the long-time behavior of weak solutions to the system 

u t + div(u ® u) - vAu = div ( - pi - L(Vd Vd) - 5{LAd - f(d)) ® d) , (1) 
divu = 0, (2) 

d t + u-Vd-5d-Vu- LAd + f(d) = 0, (3) 

describing the evolutionary behavior of nematic liquid crystal flows (we refer to the 
monographs [51 E] for a detailed presentation of the physical foundations of continuum 
theories of liquid crystals). Actually, system QTJ-flH]) can be seen as a simplification of 
the original Ericksen-Leslie model [7J [13], that still keeps a good level of compliance 
with experimental results. The model couples the Navier-Stokes equation ([1]) for the 
macroscopic velocity u (p denoting as usual the pressure), with the incompressibility 
condition (J2J) and with the equation (^ ruling the behavior of the local orientation 
vector d of the liquid crystal. Here, the function / represents the gradient w.r.t. d of 
the configuration energy F of the crystal. We choose F to be a double well potential 
having minima for \d\ — 1 and growing at infinity at most as a fourth order polynomial. 
This provides a standard relaxation of the physical constraint |d| = 1, which is very 
difficult to treat mathematically. 

In this paper, the system is complemented with the homogeneous Dirichlet 
boundary condition for u, the no-flux condition for d, and with initial conditions. 
It is settled in a smooth bounded domain Q C M. d for d = 2 or d = 3. No restriction is 
assumed on the viscosity coefficient v. 

Regarding the parameter 5, we will take 5 > 0, with the case 5 > denoting 
the presence of a stretching effect on the molecules of the crystal. Some of our results, 
however, hold only for 5 = 0. Actually, the situation 5 > is more difficult to be 
treated mathematically since the term 5d ■ Vu prevents from using maximum principle 
arguments in ()3]). For this reason, even if the initial datum do satisfies the (relaxed) 
physical constraint |do| < 1 almost everywhere, the same may not be true for d(t), for 
positive times, if 5 > 0. 

A mathematical analysis of system (JT])- (J3J) has been first addressed in the papers 
[H] and [15] (in this second work, an even more general model is taken into account). 
There, the authors consider the case 5 = and prove existence of a unique classical 
solution for d = 2, and also in dimension d = 3 under the additional assumption that 
the viscosity v is sufficiently large. These results have been extended to the case 5 > 
in the paper [TS]. Finally, the restriction on the viscosity has been recently dropped 
in [2], where weak solutions are considered and a global existence result for the 3D 
system (JTJ- (J3J) is proved in that regularity frame. Of course, uniqueness is not known 
to hold in that regularity setting. A similar result is essentially contained also in the 
recent paper (3] , where analogous estimates are derived but no formal statement of an 
existence result is provided. 

The Dirichlet boundary condition for u and either a nonhomogeneous Dirichlet 
or the no-flux boundary condition for d are treated there. Moreover, let us quote the 
recent paper [9], where these results have been extended to a more general system 
CO)-©, where also temperature effects are taken into account. We note, however, that 



the results of [5] require different boundary conditions for u (namely, the so-called 
complete slip conditions). 

The long-time behavior of system (DO)-© has been analyzed in the recent work 
[20] . still considering the case d = 2 or the case d = 3 with the large viscosity v, and 
periodic boundary conditions. More precisely, in [20] the authors show existence of a 
nonempty w-limit set for any strong bounded solution emanating from smooth initial 
data. Moreover, by using the Simon-Lojasiewicz inequality, they prove that, for the 
nonlinearity / = (\d\ 2 — l)d, this cu-limit set contains only one point. 

Stability and asymptotic stability properties of this model (actually, with even 
more complete stretching terms) have also been studied in [3], where the long-time 
behavior of solutions is analyzed in the case of periodic boundary conditions. More 
precisely, the authors prove, by means of formal estimates, that weak solutions become 
eventually smoother for large times, which suffices to have existence of non-empty 
w-limit sets. 

Finally, in the recent contribution [12], the existence of a smooth global attractor 
of finite fractal dimension is obtained in two dimensions of space. 

Our aim in this paper is to extend the results of [3J ED] in the following direc- 
tions: 

(i) we address the case d = 3 without the large viscosity assumption considering weak 
solutions; 

(ii) we consider more general C 1 functions /; 

(Hi) we use different boundary conditions and weaker initial data; 

(iv) we discuss convergence, as t tends to oo, of strong solutions in some particular 
situations. 

To get (i), we prove convergence of weak solutions, and, in some situations, we get 
strong convergence using the fact that weak solutions to the system become eventually 
smoother for times t larger than some T. This property is well-known for the (uncou- 
pled) three-dimensional N-S system, and we find conditions under which it holds also 
for the coupled system ([I])-©. Note that this result is still true for periodic boundary 
conditions, and so it improves the study done in [20]. In turn, this property (cf. (1731) - 
(175]) below) enables us to obtain properties sufficient to characterize the w-limit set. 
Assuming that / is analytic, we apply the generalized Lojasiewicz theorem to get 
convergence of the variable d. 

To address question (ii), in particular, to remove the analyticity condition, we 
make the basic observation that the set of global minimizers of the configuration energy 
of the crystal coincides with the set of constant unit vectors of M. d . Then, it is easy to 
prove that any global minimizer d satisfies the so-called normal hyperbolicity condition. 
Based on this fact, we can prove that, if the w-limit set contains a global minimizer, 
then it coincides with it (i.e., it does not contain any other point). We can also give 
two precise conditions ensuring the fact that the w-limit set contains global minimizers, 
which, unfortunately, require 6 = 0. Namely, this happens when either the diffusion 
coefficient L is large enough, or when the initial energy is very small compared with 
L (in particular, the initial datum d is already close enough to the set of global 
minimizers in a suitable norm). 



The paper is organized as follows. In the next section, we present our assump- 
tions, state the main results, and, for the reader's convenience, we briefly sketch the 
basic estimates at the core of the existence proof. The proofs of the new results on the 
long-time behavior are given in Section |3j 

Acknowledgment. The authors are grateful to Professor Eduard Feireisl for 
illuminating discussions regarding the results proved in this paper. 

2 Main results 

We let O be a smooth, bounded, and connected domain in M d , d £ {2,3}, with the 
boundary I\ For simplicity, we also assume \Q\ = 1. We set H := L 2 (Q), H := L 2 (Q) d , 
and denote by (•, •) the scalar product both in H and in H and by || • || the related 
norms. Next, we set V := H 1 ^), V := H 1 ^ and V := Hffl) d . The duality 
between V and V, as well as those between V' and V and between V' Q and Vo, will be 
indicated by (-, •). Identifying H with H' through the scalar product of H, it is then 
well known that V C H C V with continuous and dense inclusions. In other words, 
(V,H,V) constitutes a Hilbert triplet (see, e.g., [IS])- Correspondingly, we also have 
the vectorial analogues (V, H, V') and (V ,H, V ). The symbol || ■ ||x will indicate 
the norm in the generic (real) Banach space X and (•, -)x will stand for the duality 
between X' and X. 

We consider / in the form 

f(d) = (^(|d| 2 ) - l)d = \d d $(\d\ 2 ) - I^D, (4) 

where 

^GC^^+ocO^+oo)), with ^(0) = 0, ^(1) = 1 and V'(l) > 0, (5) 
is an increasing function, and the convex function ip is defined by 

f = ij, ?(1) = 1. (6) 

We also assume that there exists a constant c^ > such that 

^'( r ) < H for a11 r G [0,+oo). (7) 

Given L > 0, we define the configuration energy of the liquid crystal flow as 

£{d):=\ f (L\Vd\ 2 + $(\d\ 2 )-\d\ 2 ). (8) 

The total energy is then given by adding to £ the "macroscopic" kinetic energy; namely, 
we set 

E(«, d) := \\\u\\ 2 + 8{d) = lf (\u\ 2 + L\Vd\ 2 + ^(|d| 2 ) - |d| 2 ) . (9) 

2 2 Jq 



Let us notice that, thanks to the above assumptions flSD-([7|), £ (d) — if and only if d 
is a (constant) unit vector (cf. Lemma [2.131 below for a simple proof). 

We will address the following system of PDE's: 

u t + div(it ® u) — uAu = divS, (10) 

§ = -pi - L(Vd Vd) - 5(LAd - /(d)) ® d, (11) 

divu = 0, (12) 

d t + u- Vd-5d- Vu-LAd + f(d) = 0, (13) 

where the coefficients i/, L, 5 satisfy v, L > and <5 > 0. Notice that, by flZj), /(d) grows 
at infinity at most as the third power of \d\. 

The system, supplemented with the boundary and initial conditions 

u = a.e. on (0,T) x T, (14) 

d n d = a.e. on (0,T) x T, (15) 

w|t=o = Wo, d\ t =o = d , a.e. in Q, (16) 

will be called Problem (P). 

We introduce a precise definition of weak solutions: 

Definition 2.1. A weak solution to Problem (P) is a couple (it, d) such that 

itGL°°(0,T;if)nL 2 (0,T;V ), (17) 

d G /^(O, T; L^ 2 {Q) d ) n L°°(0, T; V) n L 2 (0, T; H 2 {Q) d ), (18) 

for all T > 0, it, d satisfy initial and boundary conditions ffTBT) . ffTBT) . the equations 
(ITT|) - (IT5|) are satisfied for a.e. t e (0, T), and 

(tit, 0) - (u®u):V<t> + v Vu:V<t>=- §:V<t>, (19) 

in in in 

holds for any test function 4> G VJ^ ' div (f2) (i.e., the subspace of W Q ' (fl) d consisting of 
divergence- free functions). 

Remark 2.2. The regularity of the test function can be justified thanks to (|17[) . 
(fT8|) and (TTSl . We have in any case (also if 5 > 0) 

it® it, VdOVd, (LAd-/(d))®dGL 2 (0,T;L 3/2 (fi)) dxd , (20) 

whence their (distributional) divergence belongs to the space L 2 (0,T; W~ 1,3 ^ 2 (^l)) d . 
Note also that the boundary condition (IT4"|) is in fact "embedded" into the weak for- 
mulation (TT9l) . 

It is known that Problem (P) admits at least one weak solution (u, d). This has been 
proved in [13] for the case 5 = and in [2] for the case 5=1 (cf. also [3J for the formal 
computations). Namely, we have 



Theorem 2.3. Let ©-© hold and let 

u eH, divu = 0, (21) 

do e V. (22) 

Then, Problem (P) possesses a global in time weak solution (u,d), satisfying, for 
a. a. t > 0, the energy inequality 

— E{u,d) + \\ - LAd + f(d)\\ 2 + u\\Vu\\ 2 <0. (23) 

CJ.6 



We point out that assumptions fl2T|) - fl22|) are equivalent to asking that the initial energy 
Eo := E(tio 5 do) is finite. 

Remark 2.4. The proof of the above theorem relies on a rather tricky approximation 
scheme and on refined compactness methods to pass to the limit. It is then worth point- 
ing out that, due to nonuniqueness, our subsequent results on the long-time behavior 
hold only for those solutions satisfying the energy inequality, in particular for the limit 
points of the approximate scheme, and not necessarily for all solutions in the regularity 
frame ffT71) - ffT8l) . Actually, there may exist "spurious" weak solutions not satisfying the 
energy inequality (T23]) which is crucial for investigating the long-time behavior. As a 
convention, in the sequel we shall restrict the terminology "weak solutions" to those 
solutions which satisfy ( 12"3"|) . Spurious solutions are thus excluded. 

As noted above, in the case 5 = a maximum principle holds for the d- 
component of any weak solution. For the reader's convenience, we recall the statement 
and the (simple) proof. 

Theorem 2.5. Let the assumptions of Theorem 12.31 hold and let 

5 = 0, and \d Q (x)\ < 1 for a. a. i6fl. 

Then any weak solution (u, d) to Problem (P) satisfies 

\d{t,x)\ < 1 for a.a. (t,x) e (0, oo) x fi. (24) 

Proof. Testing equation (JT3]) by d one obtains 

~\d\ 2 + u ■ V\d\ 2 - ^A|d| 2 + L| W| 2 + (>(|d| 2 ) - 1) \d\ 2 = 0. (25) 

Then, we notice that, by ^, 

(V>(r) - l)r > Vr > 1. (26) 

Thus, (125]) represents a parabolic equation for |d| 2 (which still satisfies the no-flux b.c). 
It is clear that the maximum principle applies, yielding fl2] 



Unfortunately, f )24|) is not known (and not expected) to hold in the case 5 > 0. 

Although the next result is essentially contained in the paper [20], for com- 
pleteness it is worth stating and proving existence of (nonempty) w-limit sets of weak 
solutions. 
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Theorem 2.6. Let the assumptions of Theorem 12.31 hold, and let (u, d) be a weak 
solution of Problem (P). Then, the co-limit set of (it, d) is nonempty. More precisely, 
we have 

lim u{t) = weakly in H, (27) 

i/"+oo 

and any diverging sequence {t n } C [0, +oo) admits a subsequence, not relabeled, such 
that 

lim d(t„) = doo weakly in V and strongly in H, (28) 

n/*-+oo 

for some d^ G V. Moreover, any sucn limit point d^ is a solution of the stationary 
problem 

- LAz + f(z) = in Q, d n z = on T. (29) 

Proof. Let {t n } C [0, +oo) be a diverging sequence. Then, the energy estimate 
implies that, at least for a (nonrelabeled) subsequence of n, 

u(t n ) — > itoo weakly in H, d(t n ) — > d^ weakly in V, (30) 

for suitable limit functions itoo and doo. Let us consider the initial and boundary value 
problem associated to (fT0|) - (fr3|) on the time interval [t n ,t n + 1] with "initial values" 
u(t n ) and d{t n ). It is clear that, setting, it n (t) := u{t + t n ) and d n (t) := d{t + t n ), 
t G [0, 1], we get a weak solution to the problem on the time interval [0, 1]. Then, (|2"3"|) 
implies that 

Vit„, -^ strongly in L 2 (0, 1; H d ), (31) 

whence, by Poincare's inequality and (1231 again, we have also 

u n -¥ strongly in L 2 (0, 1; V ) and weakly star in L°°(0, 1; H). (32) 

Moreover, we have 

d n ^d weakly star in L°°(0, 1; V) H L 2 (0, 1; /f 2 (ft) d ), (33) 

for a suitable limit function d. The growth condition ([7]) and a comparison argument 
in (fl3|) then entail 

d n>t -> d t weakly in L 2 (0, 1; L 3/2 (tt) d ). (34) 

Hence, by the Aubin-Lions lemma, we obtain 

d„ ->■ d strongly in L 2 (0, 1; V). (35) 

To proceed, we take G V^ ' div (f2) and test (fTU|) by 0. Noting that, by f TTTT) . 

[§ n :V<t>=-L f (Vd„ Vd n ) :V</>-6 [ ((LAd n - f{d n )) ® d„) : V0, (36) 
and recalling (1521 , (133|) , and Remark 12.21 we arrive at 

W U ^WmO,l;W-^ /2 (n)) ^ C ' ( 37 ) 



where W di ^ (f2) denotes the dual space to W ' div (£l) and c denotes a positive con- 
stant independent of n. 

Thus, from ([3"2"j) . ( 1371) . and the Aubin-Lions lemma, we obtain that 

u n ^0 strongly in C°([0,1]; V£), (38) 

so that, in particular, u^ = 0, and (13"2~1) . (j3"8"|) imply (1271) . On the other hand, by the 
energy estimate, we obtain 

- LAd n + f(d n ) ->■ strongly in L 2 (0, 1; if), (39) 

whereas, by (j32j)-(|55|), 

u n • Vd n - 5d n ■ Vu n ->■ weakly in L 2 (0, 1; L 3/2 (fi) d ). (40) 

Thus, comparing terms in ( JT3"|) . we also have that d n4 — >■ in a suitable way. This 
entails that d is constant in time and, therefore, it coincides with doo for all times in 
[0, 1]. Moreover, taking the limit in ( fl3|) . we obtain that d^ is a solution to (129~|) . as 
desired. This completes the proof. I 

We now present the main results of this paper, which characterize the w-limit set of 
our system as a singleton under a number of different conditions. 



Theorem 2.7. Let the assumptions of Theorem 12.31 hold, and, in addition, let ip be 
analytic. Then the oj-limit set of the component d of any weak solution consists of a 
single point, and we have 

lim d(t) = doo strongly in H (41) 

t/"+oo 

for the whole trajectory d, where d^ is a solution to (T291 • 

Remark 2.8. As usually, when applying the Lojasiewicz inequality, we can also get 
the rate of convergence of the form 

||d(t)-doo||H <C(l + t)-^s, 

where 6 is the Lojasiewicz exponent, and C is a suitably chosen constant depending 
on the initial energy and on the limit function. 

In particular situations, we can also prove a stronger convergence result: 

Theorem 2.9. Under the hypotheses of Theorem 12 .7\ let, in addition, 

5 = 0. (42) 

Then, 

u(t) -)■ strongly in V, (43) 

lim d(t) = doo strongly in H 2 {Vt) d . (44) 



Remark 2.10. The same result was proved in [20] for periodic boundary conditions 
for u and d, large viscosity coefficient, and smooth initial data. Actually, it is easy 
to check that our argument holds true also in the case of periodic B.C., when S > 0. 
Hence, the same result of [20] holds without the requirement of large viscosity, and 
for initial data as in (12T]) . f l22|) . On the other hand, in the case S > with boundary 
conditions f[H]) -f lTo^ . it does not seem possible to repeat the strong estimates required 
for the proof of (T4l?j) -( 14"4"|) (some additionally boundary terms appear, which is not 
clear how to control). Hence, extending the statement of Theorem !2.9l to this situation 
remains an open question. 

As in [2"0"| Thm. 1.2] the proofs rely on a suitable version of the Simon-Lojasiewicz 
inequality, proved in [Tf] Thm. 6]. For the reader's convenience, we report here the 
statement of a particular case of the (more general) result of [UJ, in a form suitable 
for our application: 

Theorem 2.11. Let the energy functional £ be given by ([8]) with ip analytic. Let 
p G V be a critical point of £. Then there exist constants 9 G (0, 1/2), A > and 
e x > such that the inequality 

\£(v) - £(p)\ 1 - e < A|| - LAv + f(v)\\ v , (45) 

holds for any v such that 

\\v — p\\v < ei • (46) 

To apply the preceding Theorem in our situation, we have to show that the 
inequality (j45p holds for v = d(t) in a small H- neighbourhood of d^: 

Lemma 2.12. Let the energy functional £ be given by (JSJ) with ip analytic. Let 
doo G V be a solution of (|29|) . Let K, P > be constants. Then there exist e > and 
A > such that (145)) holds for any v such that 

\\v\\ v <K, \\v - d^Wn < e, and \£(v) - £(d OD )\ < P . (47) 

Proof. We argue by contradiction. Assume that there is a sequence v n such that 

\\v n \\v<K, Vn-^docinH, \£{v n ) - £(d 00 )\ < P 



and 
Then 



\£(v n )-£(d 00 )\ 1 - 9 >n\\-LAv n + f(v n )\\ v >, n = 1,2,3,... (48) 

f(v n ) -> /(doo) in V, and Av n ->• Ad^ in V. 



This implies that 

Vv n — ¥ VcZqo in H, and, consequently, v n —¥ d^ in V. 

Hence, at least for n sufficiently large, fT4B]) holds for v = v n , p = d^. Consequently, 
also (T4"5"l) is valid. This contradicts 



In the case that / does not satisfy the analyticity condition, we can show that 
the w-limit set is a singleton only in particular situations. For this purpose, we first 
state a simple property: 

Lemma 2.13. Let (151) -([71) hold. Then, d is a global minimizer of £ if and only if d is 
a constant unit vector. 

PROOF. Thanks to (J5j)-(|7|) the function r i— > ip(r) — r has a minimum at r = 1; 
moreover, ^>{1) — 1 = 0. Thus, £{d) is always nonnegative and £{d) = if and only if 
Vd = a.e. in Q and \d\ = 1 a.e. in Q, whence the claim follows immediately. I 

Our next result is of conditional type and states that, if the cu-limit set of d(t) contains 
at least one global minimizer d of the free energy, then it has to coincide with the set 
{d}. This is a consequence of the facts that the set of global minimizers of the free 
energy is a (d — l)-dimensional smooth manifold and, on the other hand, the kernel of 
the linearized operator z \-¥ —Az + ddf(d)z is also a (d — l)-dimensional manifold. 
In other words, the so-called normal hyperbolicity condition is satisfied at d, which 
implies convergence of the whole trajectory to d. 



Theorem 2.14. Let the assumptions of Theorem 12.31 hold and let us assume that 
there exist a constant unit vector d e S d_1 and a diverging sequence {t n } such that 

lim d(t n ) = d weakly in V. (49) 

i„/"+oo 

Then, u — lim d = {d} and the whole trajectory d(t) converges to d strongly in H as 
t /• oo. If, in addition, (02) holds, then d(t) -> d in H 2 (Vt) d . 

Remark 2.15. Let us note that the same convergence result for d in H 2 (Q) d holds 
true in case 5 > with periodic boundary conditions for u and d. 

The next results only hold in the case 5 = 0. Actually, their proofs rely on the maximum 
principle proved in Theorem 12.51 In this setting, convergence to a single equilibrium 
takes place if either the diffusion coefficient L in fflBT) is large enough, or the "initial 
energy" E := E(it ,d ) (cf. (J2T|) - (I22|) ) is small enough (in other words, if the initial 
datum d is sufficiently close to the set of global minimizers). Indeed, we can prove 
the following two results: 



Theorem 2.16. Let the assumptions of Theorem l2.5l .hoJd and, in particular, let 5 = 0. 
Assume that L in ffT3l) satisfies 



L > Cq, where Cq is the best constant in the Poincare-Wirtinger inequality. (50) 

Then, the u-limit set of any weak solution starting from (uq, do) consists of a single 
point (0, doc). 



Theorem 2.17. Let the assumptions of Theorem 12. 51 hold and, in particular, let S = 0. 
Assume that there exist k > and a > 1 such that 

$(r) - r > «(1 - r) a Vr-G [0,1]. (51) 

Then, there exists e > such that, if (uq, do) satisfy Eo < e, then, the co-limit set of 
any weak solution starting from (uq, do) consists of a single point (0, doo)- 

10 



3 Proofs 

All proofs will be presented in the case d = 3, the case d = 2 being clearly simpler. 

3.1 Proof of Theorem 12771 

Energy estimate. We test (TTOj) by u and ( flBl by — LAd + /(d). Performing stan- 
dard computations and using, in particular, the incompressibility constraint (TT2"j) . we 
readily obtain the energy inequality (I2"3"j) . In particular, we get that the function 
t i— ?• E(u(£), eZ(£)) is nonincreasing, whence it tends to some (finite) value K^. More- 
over, thanks to (127)1 - (BB JI . we get 

r+oo 

E OD -E = - / V(s) ds < 0, (52) 

Jo 

where T> denotes the sum of the dissipative terms, namely 

V:=\\- LAd+ f(d)\\ 2 + is\\Vu\\ 2 . (53) 

We deduce from the energy inequality (|23l) and ©-(JZD that 

uGL°°(0,oo;/f)nL 2 (0,oo; V ), (54) 

deL°°(0,oo;V), (55) 

-LAd + f(d)eL 2 (0,oo;H). (56) 

Relations (1M|) . (1351) imply 

n- Vd-5d- VuGL 2 (0,oo;L 3/2 (O) d ) (57) 

which, together with (IB^l) . yield (cf. (|T3T) ) 

diGL^O,^;^ 3 / 2 ^)^. (58) 

Application of the Lojasiewicz inequality. Our aim is to show that there exists 
T > such that 

d t eL 1 (T,oo;L 3/2 (n) d ), (59) 

which implies convergence of d in L 3 / 2 (f2) d . The pre-compactness of the trajectory in 
H then concludes the proof of Theorem 12.71 

To this end, we first realize that there exists a constant C such that 

||n(t)|| 2 < C||Vu(t)||^ fora.a. £ > 0, (60) 

where 9 e (0, |) is the same as in ( 1451) . Indeed, if ||Vn|| < 1, then (16D1 follows by the 
Poincare inequality, and if ||Vlt|| > 1, the interpolation between Vo and Vq together 
with the boundedness of u in V' gives the same estimate. 

Now, let doo be an element of the w-limit set of d. Then, integrating (|23j) from 
to +oo, we infer that 

PgL 1 (0,oo) i (61) 

11 



and, from Lemma [2.121 and (1601) . we get 

1 



i; 



V(s) ds = S(d(t)) - Sid^) + -\\u(t)\\ 2 < CV(t)W=V, (62) 



for all t > such that ( I47p holds. Denoting by A4. this set, we obtain that "D 1 / 2 G 
L 1 (A^), see pUl Lemma 7.1]. Then also 

ueL\M;V Q ), -LAd+ f(d) e L X {M]H), (63) 

and so, taking into account the growth of /, we get 

u-Vd-8d-Vuel}{M, Ll(tt) d ), 

which implies (cf. ( TT3"j) ) 

d t eL l (M,L*(n) d ). (64) 

This fact, combined with the pre-compactness of the trajectory of d in H and a 
simple contradiction argument (see [10]), yields the existence of a large T such that 
(cf. Lemma [2TT2J 1 

||d(t)-doo||<e, Vt>T. 

In other words, the solution d remains in the e— neighbourhood of d^ in the space 
H for t > T, and Lemma [2.121 applies to d(t) in the whole interval (T, oo). In other 
words, M. D (T, oo), and (|59|) holds. This fact, together with pre-compactness of the 
trajectory yields 

d(t) — > doo strongly in H . (65) 

Theorem 12.71 has been proved. 

3.2 Proof of Theorem [2T91 

The proof follows the lines of the argument developed in [31 Sec. 3.2]. In particular, 
we can derive the differential inequality for the dissipative term T> defined in ( 153]) (cf. 
[21 formula (8)], which is still valid with our boundary conditions in case S = or with 
periodic boundary conditions in case 5 > 0): 

^P<a(P 3 + l), (66) 

where the computable constant C* depends on the parameters of the problem and on 
the "initial energy" E , but is independent of time. We point out that the formal 
computations in the proof of (1661) are valid for the classical solutions, but the result 
can be justified by a proper approximation. 

Next, we show that there is T\ > such that T> G L°°(Ti, oo). To prove this, we 
consider the differential inequality 

y'<a(y 3 + l), y(t ) = l. (67) 
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Then, there exist (a small) r (independent of to) an d (a large) K > such that the 
solution y satisfies 

IMIc°([to,to+r]) < K. (68) 

On the other hand, according to (ISTj) . we have 



Thus, for any e > there exists T > such that 



r+oo 

lim / Vis) ds = 0. (69) 



V(s) ds < e. (70) 



r 



Choosing e = r/2 and T correspondingly, we obtain that, for all t > T, there exists 
to ^ [t, t + r/2] such that 

2 f t+T/2 2e 

V(to) < ~ / v (s) ds < — = 1. (71) 



Comparing solutions of f l66|) and fl67|) and recalling the choice of to, we get from fl68|) 
that 

\V{s)\<K VsG [T + r/2,+oo). (72) 



Setting Tx := T + r/2, we deduce from ( 16"B"|) that ^-P is bounded on (T l5 oo), which 
together with (16TT) yields 

X>(t) ->• as t -> oo. (73) 

This implies (using the Poincare inequality) that, 

u(t) -)> 0, strongly in V. (74) 

Taking into account the growth of /, we get from fT72|) 

MllL°°(Ti,oo;^ 2 (f7) d ) < -^- (75) 

To show that d converges to a single point d^, we make again use of the Lojasiewicz 
inequality ( 145|) . The same argument applies this time to the strong solution and time 
t>T\. This gives 

d t e I}(T, oo; H) for some T > T x . (76) 

It follows that 

d(t) — y doo strongly in H. (77) 

Moreover, by (17511 and the growth conditions on /, 



f(d(t)) -> /(doo) strongly in//, (78) 

whence, using ( TT3|) and ( f74|) . 

Ad(t) -> Adoo strongly in H, (79) 

which concludes the proof. 
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3.3 Proof of Theorem [2341 

In this section, we show that the energy functional £ satisfies the Lojasiewicz inequality 
fj45|) with the exponent |. Then, arguing as in the proof of Theorem 12 .71 we obtain the 
strong convergence of d in H . Moreover, if (j4~2j) holds (or we have periodic boundary 
conditions, cf. Remark 12.101) . we have the strong convergence in H 2 (Q) d (cf. the proof 
of Theorem 12. 9p . 

Let us consider the linearized problem associated to (12T?|) at the element d of the 
w-limit set, i.e., 

£(d)z := -LAz + ^(\d\ 2 )z + 2ij'(\d\ 2 )(d®d)z- z = 0, d n z = 0onT. (80) 

Let d be a global minimizer of £, i.e., a constant unit vector (by Lemma l2.13p . We 
aim to apply the result proved by Simon and reported in jH Cor. 3.12]. To this end, 
we introduce the following notation: 

U is a "^-neighbourhood of d G V, 

Vo is the kernel of C(d), 
S = {de V; £'{d) = 0}, 
S = {heU; £'(d + h) eV }. 
In our situation, jU Cor. 3.12] reads as follows: 

Lemma 3.1. Let d G Sq and assume the following hypotheses: 

(i) The kernel Vo of the linearization C(d) is a complemented subspace of V, 
i.e., there exists a projection P G £>( V) such that Vo = Rg-P- 

(ii) There exists a neighbourhood U of d in V such that £' G C 1 (t/, V). More- 
over, the range of C(d) coincides with V[, the space of the elements of V' belonging 
to the the kernel of the adjoint projection P' G B(V'). 

(in) (S — d) Pi S is a neighbourhood of in the critical manifold S. 
Then £ satisfies the Lojasiewicz inequality near d with the exponent 9 = \. 

To verify the assumption (i), we test (18"U|) by z and use the condition -0(1) = 1 to 
obtain 

L||Vz|| 2 + 2 f ij\l)\d-z\ 2 = 0. (81) 

Jn 

Hence, taking into account the last condition in ([5]), we get 

Vz = and d ■ z = a.e. in fi. (82) 

Consequently, any solution z to (IHUj) . i.e., any element of the kernel, is a constant 
vector orthogonal to d (conversely, it is apparent that any such vector is a solution 
to (I80p ). Thus, the kernel of the linearized operator C(d) is a (d — 1) -dimensional 
plane orthogonal to d and containing the origin, which trivially permits to define the 
projection P. 

The first condition in (ii) is obvious since / is C 1 and, by hypotheses, has at most 
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cubic growth. To verify the second condition, we observe that Vi is the subspace of 
V' consisting of the elements that are orthogonal (w.r.t. the duality between V' and 
V) to the plane Vo- Then, computing (£(d)z, v ) for generic z <E V and v G Vo, we 
obtain (cf. (JHQJ and recall that ip(l) = 1), using (l8"2j). 



(C(d)z, vo) = 2 / ^'(|d| 2 )(d -z)(d- v ) = 0, (83) 

Jn 

the last equality following from the fact that v Q _L d. Thus, C(d)V C Vi- To show the 
converse inclusion, we choose £ G V^ and prove that there exists at least one z G V 
such that 

C(d)z = -LAz + 2ip'(l)(d®d)z = C, d n z = 0onT. (84) 

This can be seen by approximation. Actually, it is clear that, for any fcsN, there is a 
solution z k to 

- LAz k + k~ 1 z k + 2tjj'(l)(d®d)z k = £, 9 n z = on T. (85) 

Testing by z k , we have 

L||V2 fc || 2 + k-'WzkW 2 + 2ij'(l)\\d-z k \\ 2 = (C, z k ) = (C, z k - Pz k ). (86) 

Indeed, C, -L -Pz/c by assumption. Using the Poincare-Wirtinger inequality it is then 
apparent that the right-hand side can be estimated. Then, standard methods permit 
to check that z k tend to a solution z to (1541) . as desired. 



The third assumption is satisfied because G So — d C S, and both So — d and S have 
the same dimension (see [H Proposition 3.6]). 

Lemma 13.1 1 then yields that the Lojasiewicz inequality (|4"5"|) holds near d with 
the exponent 0=1/2. If 5 = or in case of periodic boundary conditions, repeating 
the computations leading to (175]) . we obtain the strong convergence of d in H, and, 
proceeding as in (T78|) . (T79|) . the strong convergence in H 2 (Q) d , which completes the 
proof of Theorem 12.141 

3.4 Proof of Theorem [2361 

Let doo be an element of the cu-limit set of the d-component of some weak solution. 
Then, by Theorem 12. 6[ d^ solves the stationary problem, which we rewrite as 

-LAd oo +(i;(\d oo \ 2 )-l)d oo = 0. (87) 

Testing (JHZD by doo - (doo)fi, where (d^ = J n d^, we get 

LUVdoof + / V(|doo| 2 )doc ■ (doo - (doo)n) = ||doo - (rfoo)n|| 2 < c^UVdooH 2 , (88) 
in 

where cq is the (best) constant in the Poincare-Wirtinger inequality. Thus, being L is 
large (precisely, we need L > Cq), the latter term can be controlled. 
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In what follows, we denote 

5(d) = \mn (89) 

A direct check (e.g., computing the Hessian matrix) shows that S is convex, and we 
notice that 9 d 5(d) = ip(\d\ 2 )d. Thus, we have 

/ V(|doo| 2 )doo • (doo - (doo)«) = (9d5(doo),doo - (^00)^)^ 

> / (S(doo) - S((doo)n)) 
./n 

= JEid^-E^ /doo) >0, (90) 

where the latter inequality follows from Jensen's inequality. 

From (J88p - (j90p . we obtain that doo is a constant vector. Thus, taking into 
account that ^ is monotone and "0(1) = 1, we readily obtain from equation ( 187)1 that 
either d^ = or |doo| = 1- Hence the set of stationary solutions is disconnected and 
consists of the isolated point and the two-dimensional manifold. Consequently, either 
Theorem 12.141 applies, or the whole trajectory tends to 0. In both cases, the cu-limit 
set is a singleton, as desired. 

3.5 Proof of Theorem 12371 

Let us first note that, by (1521) and (15TJ) . we have 

e > E > E^ = ~ / (L|Vdoo| 2 + ^(|doo| 2 ) - |doo| 2 )- (91) 

^ Jn 

Rewriting the stationary problem ([57]) and testing it by d^, we obtain 

£||Vdoo|| 2 + / (^(IdooHldool 2 - (dool 2 ) = 0. (92) 

Jn 

Dividing (I9"2~j) by 2 and subtracting the result from (19~TJ) . we obtain 

\ /* (^(Irf^l 2 ) - ^(|rfoo| 2 )|rfoo| 2 ) < e. (93) 

2 Jn 

On the other hand, thanks to dSJ, © and (15T|) . 

i / (^(Id^l^-^dd^l^ldool 2 ) >\f (^(Id^l^-ldool 2 ) > £ / ll-ld^l 2 ^, (94) 

1 Jn z Jn *> Jn 

where also the maximum principle (|24p has been used. Thus, 



'2e\iA 



|1- |*o| 2 |L,, m < 111 - |*o| 2 |Lw m < - • (95) 
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To proceed, we notice that, by standard elliptic regularity results applied to (I2"9~j) . there 
exists a constant K > such that, for any solution d of (|2"9"|) it holds 

IfflUn) < *o- (96) 

Consequently, for some K > depending on K , we have 

Ili-H'lUo)^^- ( 97 ) 

In particular, d^ satisfies (19~T|) . Thus, by the Gagliardo-Nirenberg interpolation in- 
equality (we refer, for simplicity, to the case d = 3, the case d = 2 is even better), 

111 -Id l 2 ll < 111 -Id l 2 ll 6/7 111 -Id l 2 ll 1/7 + 111 -Id l 2 ll 

where rj > is a small constant to be chosen later and the last inequality holds provided 
that e is small enough. 

To conclude the proof, we set 

a : = ^(Idool 2 ) - 1 (99) 

and notice that a < because of ©-©• Moreover, doo can be interpreted as a solution 
of the linear elliptic system 

- LAd^ + ad^ = in Q, d^d^ = on T. (100) 

Testing the above equation by 1, we obtain 

adoo = 0. (101) 

Then, multiplying fllOOp by d^ — (d^n and using fllOll) . we infer 

£||Vdoo|| 2 = - / ajdoo - (doo)^ 2 - (doo)^ • / a(d OQ - (d^n) 
Jn Jn 



<^j||a|U»(n)||Vd 00 || +|(doo)n| a. (102) 

Jn 

The latter term is nonpositive, while the first term on the right-hand side can be 
controlled provided that r\ is small enough. Indeed, 

NU»(n) = ||^(l d oo| 2 ) - l|| L oc (n) = \\^(\doo\ 2 ) - ^(l)|| L oc (n) 

< c^HlcZool 2 -l|| Loo(n) <c^r] } (103) 

thanks to (19"8|) and ([7j). If L > c^c^rj, we see, in the same way as above, that d^ is 
a constant unit vector. Finally, we take e such that (1951) holds, and the proof follows 
again by applying Theorem 12.141 
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